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ABSTRACT Translation 



The invention concerns a method for calculating the relative volumetric flow-rates of at 
least one of the phases of a multiphase effluent flowing in a well, said method comprising 
a first step of acquiring local volumetric fractions and/or velocities of said phases across a 
section of the well at a certain depth. The method further comprises: - corrections of said 
local volumetric fraction and/or velocity measurements in order to make them consistent 
with each other and/or with the effluent flow conditions; - selection of a suitable flow 
model mathematically representing the effluent flow; - interpolation of said local 
volumetric fraction measurements and/or said local velocity measurements by the 
selected flow model in order to obtain a volumetric fraction profile and/or a velocity 
profile for at least one phase of the effluent across said section of the well at said depth.- 
calculation of the relative volumetric flow-rates of said at least one phase by integration 
of said volumetric fraction and/or velocity profiles over said section of the well at said 
depth. 
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© Precede d'analyse d'un ecoulement diphasique dans un puits d'hydrocarbures. 

© L' invention a pour objet un procede d'analyse 
de I'ecoulement dans un puits incline, recoup- 
ment comprenant du petrol e et de I'eau. 

Le procede selon ('invention est caracterise 
par le fait que, pour une profondeur ou Ton 
dispose d'au moins une mesure d'un parametre 
d'ecoulement, on eta b I it a partir de cette 
mesure, connaissant I'inclinaison 6 et le dia- 
metre D du puits a la dite profondeur, un profil 
de variation de ce parametre le long d'un axe z'z 
defini par le diametre de la section transversale 
du puits situe dans un plan vertical, en utilisant 
un modele prevoyant une couche de melange 
eau-petrole turbulente. 
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[.'invention concerne de facon generate I'analyse d'ecoulements diphasiques eau- petrole dans tes puits 
d'hydrocarbures. L'invention concerne plus particulierementun precede permettantde determiner, a partirde 
mesures, faites a des profondeurs donnees, de parametres tels que ia proportion statique d'eau, un profit de 
variation dudit parametre suivant la section transversale du puits a la profondeur consideree. L'invention 
5 concerne aussi la determination du debit de petrole a partir de tels prof ils. 

II est habitual dans la production d'hydrocarbures que le fluide produit par le puits soit un melange de deux 
phases, petrole et eau. Le petrole constitue le plus frequemment la phase d is perse e, formee de bulles, et I'eau 
la phase continue. 

Les techniques classiques pour I'analyse des ecoulements dans les puits d'hydrocarbures reposent sur 
10 la mesure du debit total Q t et de la proportion statique moyenne d'eau H w (en anglais "water holdup"). On utilise 
habituellement un debitmetre a he! ice, sensible a la vitesse, pour mesurer le debit total Q t , et un appareil a 
pression differentielle appele gradiomanometre, sensible a la densite, pour la proportion statique d'eau 
On determine le debit de petrole Qo a partir de ces mesures par la formule classique 

Qo = (1 - H w ) (Q t + V s .rVS) 
15 ou V s est la vitesse de glissement entre phases et S la section du puits. 

On fait ici I'hypothese que les valeurs de vitesse et de proportion statique resultant des mesures precitees, 
qui sont des valeurs moyennes sur la section d'ecoulement, sont bien representatives de I'ecoulement, en d'au- 
tres termes, que les variations de la vitesse et de la proportion statique d'eau sur une section du puits sont 
negligeables, ou en tout cas ont globalement une influence negligeable sur le resultat. 
20 Cette hypothese est justhf iee pour un puits vertical. Mais lorsque le puits est incline, la difference de densite 

entre les phases tend a provoquer une segregation des phases : I'eau, phase lourde, predomine dans le bas 
de la section d'ecoulement, tandis que le petrole, plus leger, predomine dans le haut de cette section. Dans 
certaines conditions, la segregation des phases est completement realisee, a savoirque le bas de la section 
d'ecoulement est uniquement occupe par la phase eau et le haut par la phase petrole (voir a ce sujet I'article 
25 SPE 1821 7 "Oil/Water Flow Structure Measurements in Inclined Pipes". 2-SOctobre 1988, de P. Vigneaux, G. 
Catala et J.P. Hulin). On observe dans ce type de situation une circulation a contresens ("backflow") dans le 
bas de la section d'ecoulement 

II est evident que dans ces conditions, les variations de la proportion d'eau et de la vitesse selon la section 
d'ecoulement etant importantes, les resultats obtenus paries mesures precitees ne fournissent pas une esti- 
30 mation f iable des valeurs moyennes recherchees. 

L'invention a pour objet un procede d'analyse de l'6coulement dans un puits incline, I'ecoulement compre- 
nant du petrole et de I'eau, lequel procede est def ini dans les revendications annexees. 

Le profil de variation obtenu selon l'invention represente une information precieuse pour l'interpr6tation 
des mesures. II peut etre etabli a partir d'une pluralite de mesures locales telles que celtes decrites dans le 
35 brevet US 4 974 446 ou la demande de brevet europeen 0 394 085, et permet alors de determiner une valeur 
moyenne par integration sur la section d'ecoulement. 

II peut egalement etre etabli, en sens inverse, a partir d'une valeur moyenne, obtenue par une mesure de 
type classique telle que decrite plus haut. 

L'invention sera bien comprise a la lecture de la description ci-apres, faite en reference aux dessins an- 
40 nexes. 

Dans les dessins : 

- la figure 1 represente un model e d'ecoulement utilise selon l'invention; 

- les figures 2Aa 2D represente nt les d if fe rents types de profits de variation de la proportion statique d'eau 
a, dans le cadre de ce modele; 

45 - les figures 3 et 4 illustrent la concordance des valeurs experimental es avec ceiles predites par le second 

modele precite; 

- la figure 5 est un diagramme decrivant I 'eta bliss em ent d'un profil de variation de a a partir de mesures 
locales de ce parametre; 

- la figure 6 est un diagramme decrivant I'etablissement d'un profil de variation de a a partir d'une mesure 
so de la valeur moyenne H w de la proportion d'eau; 

- la figure 7 represente schematiquement une section d'ecoulement du point de vue des vitesses; 

- la figure 8 est un diagramme decrivant I'etablissement d'un profil de vitesse a partir d'une mesure de la 
valeur moyenne de la vitesse; 

- la figure 9 est un diagramme decrivant I'etablissement d'un profil de vitesse et du debit de petrole a partir 
55 de mesures locales de vitesse; 

- la figure 10 permet de comparer des valeurs de vitesse experimentales avec les valeurs predites par le 
modele. 

La figure 1 represente une portion d'un puits de production de petrole de diametre D. Dans ce puits circule 

2 




25 



vers la surface un fluide forme d'un melange de petrole et d'eau de formation. On supposera dans ce qui suit 
que le fluide est diphasique et que ie gaz present est en solution dans le fluide. Le cas le plus frequent est 
que le petrole constitue la phase disperses et I'eau la phase continue et c'est a cette situation que I'on se re- 
ferera, etantentendu qu'il ne s'agit pas d'une limitation de I'invention. 

5 On def in it un systeme de coordonnees comprenant un axe x'x situ e sur I'axe du puits, et des axes y'y (non 

represente) et z'z situes dans un plan perpend iculaire a I'axe x'x, I'axe z'z etant dans un plan vertical. L'axe 
x'x est incline d'un angle 9 par rapport a la verticale. 

Les deux phases ont des densites differentes, de I'ordre de 0,8 pour le petrole, entre 1 et 1,1 pour I'eau, 
habituellement salee. II en resulte une tendance a la segregation des phases : le petrole tend a predominer 

10 dans le haut de la section d'ecoulement, et I'eau a predominer dans le bas de la section. 

La figure 1 illustre le model e utilise selon I'invention pour I' etude de cet ecoulement. Le model e suppose 
I' existence d'une couche turbulente de melange eau-petrole, dont I'extension selon I'axe z'z est notee X. La 
figure 1 montre en trait epais un prof il de variation de la proportion statique d'eau a (en anglais "water holdup"') 
selon I'axe z'z : dans I'exemple represente, la proportion d'eau vaut 1 dans le bas de la section, 0 dans le haut, 

15 et varie dans la couche de melange avec un gradient constant. Le prof il de la figure 1 correspond a un modele 
simple, entierement def ini par la dimension X et par la position Zc du centre de la couche de melange. La pro- 
portion d'eau a est alors don nee par : 

° a m 1 pour z < z c - ^ (la) 

Z 2 

a - J + pour z c -2<z< z c + 2 

a m 0 pour z > z c + | (lc) 
Si Ton appelle a m la moyenne de la proportion d'eau a suivant I'axe z'z, I'expression (1 b) ci-dessus devient : 

30 z 

a = a m -f (2) 

L'exemple de la figure 1 montre une couche de melange entierement contenue a Tinterieur de la section 
d'ecoulement Mais le modele prevoit plusieurs situations selon les valeurs de X et z^, qui sont illustrees a la 
figure 2. 

35 La figure 2A montre un profit du meme type que celui de la figure 1, a savoir que la couche de melange 

est entierement contenue dans la section d'ecoulement. Le fluide dans le bas de la section est uniquement de 
i'eau, et seule la phase petrole est presents dans le haut de la stction. Cela correspond a une dimension X 
plus petite que le diametre D, et a une position Zc proche de I'axe x'x. 

La figure 2B montre un prof il ou le bord inferieur de la couche de melange est a I'interieur de la section 

40 d'ecoulement mais ou le bord superieur est a I'exterieur. Dans ce cas, seute la phase eau est presente dans 
le bas de la section tandis que le fluide occupant le haut de la section est un melange des phases a faible 
proportion d'eau. Ceci correspond a une dimension X plus petite que le diametre D eta une position positive 
et fortement ecartee de I'axe x'x. 

La figure 2C correspond au cas ou la dimension X de la couche de melange est superieure au diametre 

45 de la section. Toute la section est occupee par un melange des phases eau et petrole. 

La figure 2D montre un prof il analogue a celui de la figure 2B mais en sens inverse. La dimension X est 
plus petite que le diametre D, et la position z<. est negative et fortement ecartee de I'axe x'x. Seule la phase 
petrole occupe le haut de la section d'ecoulement, alors que le bas de la section est occupe par un melange 
a forte proportion d'eau. 

50 En utilisant une analyse de dimension, et en negligeant la viscosite, on aboutit a la relation ci-apres pour 

la dimension X: 

X = — R— (3) 

ou K est une constante pour un diametre donne (egale a 1 ,06 pour D = 20 cm et 1 ,4 pour D = 10 cm). 
55 Pour etablir le prof il de variation a(z), il faut en outre determiner le parametre z c ou, ce qui revient au meme, 

la moyenne a m de la proportion d'eau selon i'axe z'z. Ceci est realise en ajustant le parametre Zc ou a m de 
facon a faire concorder les valeurs calculees suivant le modele avec des valeurs mesurees. On utilise pour 
cela des mesures de la proportion statique d'eau a a la profondeur consideree, qui peuvent etre soit une ou 
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plusieurs mesures de la valeur locale de a, soit une mesure de la valeur moyenne de a sur la section d'ecou- 
lement 

Le modele geometrique expose ci-dessus presente dans la plupart des cas une concordance raisonnable, 
compte tenu de sa simplicity avec les donnees experimental disponibles, presentees dans Tarticle precite 
5 de Vigneaux et al. Un modele plus precis a ete etabli, selon une approche fondee sur ies lois de la diffusion 
turbulente. 

Un element du melange, d'epaisseur dz et de densite p(z), est le resultat de deux flux opposes : Tun du a 
la force d'Archimede liee a la difference des densites, qui tend a faire monter la phase petroie, plus legere, 
Tautre du a la diffusion turbulente, qui tend a homogeneiser le melange. 
10 Cela s'exprime par ('equation 

K fe + s(1 - a) = 0 (4) 

ou K est un coefficient de diffusion turbulente et s est la vitesse de sedimentation. En supposant s constant, 
on aboutit a ['expression 
15 a(z) = 1 - C.e^ (5) 

II est avantageux d'exprimer a(z) en fonction de la moyenne a m de a(z) selon I'axe z'z. Considerons le 
cas d'un prof il du type represents a la figure 2C, ou le bord superieur de la couche turbulente est exterieur a 
la conduite. II n'existe alors aucune valeur de z comprise entre -D/2 et +D/2 pour laquelle a = 0. Par integration 
sur Tintervalle {-D/2, +D/2}, on trouve 

a(Z) = 1 - P(1 s - p am W (6) 



avec 



Z = (grandeur sans dimension) 



25 _ SD 

2K 

En confrontant le gradient calcule a partir de I'expression (6) au gradient constant du premier modele, on 
obtient pour P I'expression 

p _ K.sin8 m 

30 2(1 - W 

K ayant la m£me valeur que dans le premier modele, ce qui donne 

a(Z) = 1 - ^^e2P2 (8) 
' 2shP x ' 

L'expression (6) ci-dessus a ete obtenue en supposant a(z) > 0 sur tout le prof il. Dans le cas d*un prof il 
35 du type represents a la figure 2D, ia proportion d'eau devient nulle a partir d'un certain point Cela correspond 
a I'existence d'une valeur z« qui remplit simultanement les conditions : 
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I < D/2 



Le prof il a(z) se compose alors d'une portion de courbe sur Tintervalle {-D/2, zj et d'une portion a = 0 
45 sur Tintervalle fe, +D/2}. On a constate que l'expression (6), bien qu'elle corresponds en toute rigueur a un 
cas different, a savoir > D/2, four nit une approximation satisfaisante pour la portion du prof il a(z) dans Tin- 
tervalle {-D/2, zj. 

On peut alors, de facon generale, presenter la definition du prof il comme suit : 
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K.sinO 



ot(Z) - 1 - 



2PZ 



2shP 



pour z < min{z , D/2} 



et, uniquement dans le cas ou z^ < D/2 



«<2) 



pour < z < D/2 
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min{Ze, D/2} signif iant celle des valeurs de et D/2 qui est la plus petite. 

Le prof ii a(z) est etabli en determinant la valeur du parametre P (ou a m ) par ajustement du modele en fonc- 
tion de mesures de la proportion d'eau a. 

Les figures 3 et 4 illustrent la bonne concordance de ce modele (courbes en trait plein) avec les valeurs 
5 experimental (pointilles) obtenues dans une installation a conduite inclinable. Les conditions sont les sui- 
vantes. 

Dans le cas de la figure 3 r I'inclinaison e est 65°, le diametre D est 200 mm, le debit total 2000 barils/jour 
(1 3,2 rr^/h). Les quatre courbes reperees A1 , A2, A3, A4 correspondent respectivement a des proportions dy- 
namiques d'eau (en anglais "watercut") de 95%, 90%, 70% et 30%. Le parametre porte en abscisse est la pro- 
10 portion statique d'eau a, le parametre en ordonnee est le rapport z/R (R rayon de la conduite). 

Dans le cas de la figure 4, la proportion dynamique d'eau est f ixee a 50%, et les courbes correspondent 
a des valeurs drfferentes pour I'inclinaison 6 : 65° pour la courbe B1, 25° pour la courbe B2, 5° pour la courbe 
B3. 

Comme indique plus haut a propos du modele a gradient constant, on peut, pour etablir ce profil, utiliser 
15 soit des mesures locales de a, soit une mesure de la valeur moyenne sur la section d'ecoulement. Les etapes 
de calcul utilisant les mesures locales sont illustrees par le diagramme-blocs de la figure 5, et le calcul a partir 
de la valeur moyenne est represents par le diagramme-blocs de la figure 6. 

L'instrument decrit dans le brevet US 4 974 446 ou la demande de brevet europeen 0 394 085 possede 
un nombre N de sondes hyperfrequence de petite dimension, angulairement reparties sur un cercle centre sur 
20 I'axe de l'instrument, de diametre D t , la position angulaire d'une sonde par rapport a une origine liee a I'axe 
z'z etant notee p,(i=1 ,N). L'instrument peut etre assujetti a rester centre par rapport au puits, ou bien concu de 
facon a prendre appui sur la partie inferieure de la paroi du puits. 

Des moyens sont associes a l'instrument pour mesurer I'inclinaison 8 du puits et les positions angulaires 
5|. Le diametre D du puits est suppose connu, ou mesure a Taide d'un diametreur. 
25 Les valeurs de D et 0, ainsi que de la constante K, connue pour un diametre D donne, sont fournies dans 

l'etape 51 du diagramme. 

L'etape 52 est I'entree des donnees de mesure fournies par l'instrument susvise, a savoir, le diametre D t , 
la mesure locale de proportion d'eau ai et les positions angulaires respectives 8| des sondes. 

L'etape 53 fournit pour chacune des sondes i la coordonnee Z| selon I'axe z'z, en fonction du diametre D t 
30 et de la position angulaire 5 ( . Le calcul depend bien entendu de la facon dont ('instrument est place par rapport 
au puits (centre, ou applique vers le bas). 

Dans I'etape 54, on determine une valeur initiale pour le parametre <x m , en calculant la moyenne des me- 
sures Oj. 

On calcule ensuite (etape 55) par la methode des moindres carres la valeur a™ qui correspond a I'ajuste- 
35 ment optimal de I'ensemble des valeurs a(zi) calculees comme indiqu6 plus haut aux valeurs de mesures aj. 

Ce resultat a m est la moyenne de la proportion d'eau a selon I'axe z'z. On en tire le profil a(z), a partir de 
quoi on peut calculer (etape 56) la proportion moyenne d'eau H w , par integration du profil a(z) sur la section 
d'ecoulement circulate A, selon ('expression 



"2 
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Des representations graph iques du profil a(z) et de la proportion moyenne d'eau H w peuvent Stre obtenues 
en sortie. 

Le diagramme de la figure 6 montre les etapes du traitement permettant d'etablir un profil a(z) a partir de 
la proportion moyenne d'eau H w . Dans l'etape initiale 61, sont fournies les valeurs du diametre D, de I'inclinai- 
so son 9, de la constante K, et la valeur de obtenue par mesure. 

On fixe dans l'etape 62 une valeur initiate (aJi pour a m , qui peut Stre prise egale a H w . 

L'etape 63 consiste dans le calcul d'une valeur moyenne (H) jt prise sur la section d'ecoulement de forme 
circulaire. Ce calcul se fait conformement a ('expression (9) ci-dessus. 

Cette valeur (H)j est comparee a dans l'etape 64. Si la difference n'est pas inferieure a une quantite e H 
55 determinee, le calcul est recommence avec une valeur (ct m ) j+1 pour a m differente, et ce jusqu'a ce que la dif- 
ference entre le resultat de l'etape 63 et H w soit inferieure a la quantite 6 H . La valeur correspondante de <x m 
per met d'etablir le profil a(z). 

On se refere ensuite a la figure 7 qui montre un profil de vitesse v(z), c'est-a-dire la variation de la vitesse 
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v suivant i'axe z'z def ini comme precedemment (axe perpendiculaire a I'axe du puits x'x et situe dans un plan 
vertical). La vitesse v considered ici prend en compte I'ensemble de I'ecoulement, c'est-a-dire les deux phases. 
C'est la vitesse du melange : 

v = a.v w + (1 - a).v G 
5 ou v w et v 0 sont respective me nt les vitesses des phases eau et petrole. 

Le prof il represents a titre d'exemple sur la figure 6 comporte une partie negative dans le bas de la section. 
Cela correspond a un ecoulement vers le bas du fluide occupant la partie inferieure de la section, tandis que 
le fluide occupant le haut de la section s'ecoule vers la surface (partie positive du prof il). 

Du point de vue theorique, le prof il de vitesse v(z) est la solution de I'equation de Navier-Stokes, dans le 
10 cas d'un gradient de densite p(z) donne par 

p(z) = pwa(z) + p 0 (1-a(z)) (10) 
Pw et p 0 etant les densites respectivement de I'eau et du petrole prod u its dans le puits. 

En falsant les hypotheses simplificatrices d'un profil a(z) correspondant a I'expression (1b) ci-dessus, et 
d'une viscosite turbulente constante, et en supposant, en ce qui concerne les conditions aux Nmites, que la 
15 vitesse au voisinage des parois obeit a un profil logarithmique, on aboutit a I'expression 




v(z) - - 6V m (Z2 - t) + X .D <sine.cose) % .2 

- §S sin*8.Z(Z* - i) (11) 

ou 

V m est la vitesse moyenne selon I'axe z'z 
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X est une constante egale a environ 5.s _1 

0 est une constante egale a environ 0,7 m 2 s~ 1 

Dans ce cas, la valeur de P, def ini a I'expression (7), est relativement grande, comprise entre 1 et 50 ou 
30 plus, par exemple. P est fonction de la proportion d'eau et une valeur relativement importante indique des gra- 
dients importants et non-uniformes. 

Le profil de vitesse v(z) peut etre etabli a partir d'une mesure de la valeur moyenne de la vitesse sur la 
section d'ecoulement, telle que la four nit un debitmetre a helice, ou bien a partir de mesures locales de la vi- 
tesse des bulles ou de la vitesse superf icielle de la phase dispersee, obtenues conformement a I'enseignement 
35 de la demande de brevet europeen precitee 0 394 085. 

Les traitements correspondants sont representes respectivement par les diagrammes des figures 8 et 9. 

Le traitement utilisant la vitesse moyenne comprend une etape initiate 81 dans laquelle sontfournies les 
valeurs des para metres D et 6 (diametre et inclinaison du puits), des constantes p et x, et de la vitesse moyenne 
V (egale au rapport du debit total Q T sur la section A), mesuree a la profondeur consideree. 
40 Le calcul iteratif qui suit permet d'ajuster la grandeur V m , moyenne de la vitesse v(z) suivant I'axe z'z, de 

maniere a faire concorder la vitesse moyenne correspondante sur la section d'ecoulement avec la valeur me- 
suree. 

L'etape 82 consiste a fixer pour V m une valeur initiate (V m ) 1( que I'on choisit egale a la valeur mesuree V. 
Dans l'etape 83, on calcule la valeur moyenne de la vitesse sur la section d'ecoulement de diametre D, 
45 selon la relation 



50 



r 



v(z){V m = V}.<D' 



4z*) 1 \dz 



(12) 



ou la notation v(z){V m = V h ) signif ie que Ton donne au parametre V m figurant dans I'expression de v(z) la valeur 
V. 

55 La valeur (V)j est comparee a la valeur mesuree V dans l'etape 84. Si la difference est superieure a une 

petite quantite determinee e v , on poursuit I'iteration en modrf iant la valeur donnee a V m (etape 85), et ce jusqua 
ce que la difference soit inferieure a la quantite s v , auquel cas on retient la valeur correspondante (V m ) k comme 
valeur de V m , et I'on etablit le profil de vitesse v(z). 
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Lediagrammede la figure 9 indiquetes etapes du calcul pourl'etablissementdu profil de vitesseen partant 
de mesures locales. 

Les mesures locales obtenues conformement a I'enseignement de la demande de brevet europeen 0 394 
085, au moyen d'un instrument equipe de N sondes hyperfrequence, permettent d'obtenir pour chaque sonde 
5 notee i (i=1 ,N) la proportion d'eau locale et la vitesse superf icielle vsj de la phase dispersee. Ces grandeurs 
sont liees a la vitesse locale v, par la relation 

(1 - ccO.v, = vs, - a,(1 - ai).(v 8 )j (13) 
ou (v 8 ), est la vitesse locale de glissement entre les phases. 

Selon le modele de Choquette, la vitesse locale de glissement est donnee par la relation 
10 vs = V„ m .e- «0 ■ «> (14) 

ou c est une constante et V„ m est de la forme ffp*, - p o ).(cos0) 1 /2. 

L'etape initiale 91 fournit les valeurs des constantes x et p mentionnees ci-dessus en reference a I'expres- 
sion (11), et des parametres D (diametre du puits) et 9. 

L'etape 92 consiste dans I'entree des donnees suivantes : le diametre D t du cercle sur lequel sont placees 
15 les sondes, et pour chaque sonde i, la position angulaire 6 h la proportion d'eau a, et la vitesse superf icielle 
vs,. 

L'etape 93 est identique a l'etape 53 decrite en reference a la figure 5, et fournit les coordonnees z, des 
sondes sur I'axe z'z. 

L'etape 94 consiste dans le calcul des vitesses locales v t a partir des valeurs mesurees precitees, sur la 
20 base des expressions (13) et (14). 

L'etape 95 consiste a calculer les vitesses locales v(z,) sur la base du modele correspondant a I'expression 
(11) ci-dessus, pour les valeurs z, obtenues a l'etape 93. On choisit comme valeur initiale de V m la moyenne 
des vitesses locales V| obtenues a l'etape 94. 

L'etape 96 consiste a obtenir par la methode des moindres carres la valeur de V m qui correspond au meil- 
25 leur ajustement de I'ensemble des valeurs v(zj calculees par le modele avec les valeurs v, resultant des me- 
sures. Une fois obtenue cette valeur, le profil de vitesse v(z) peut etre etabli. 

On peut evaluer les debits totaux de petrole Q Q et d'eau Q w sur la base des profits a(z) et v(z), par integration 
sur la section d'ecoulement A des debits locaux respectifs q 0 et q w donnes par les expressions 

q c = (1 - a).v + a(1 - a).v 8 
30 q w = a.v - a(1 - a).v s 

Les vitesses moyen nes correspondantes, V Q pour le petrole et V w pour I'eau, peuventetre calculees a partir 
des debits Qo et Q w : 

V = °o v = _9w_ 

° A(1 - H w ) w A.Hw 

35 Connaissant les vitesses moyennes, on peut ensuite determiner une vitesse moyenne de glissement Vs 

par: 

Vs = V 0 - V w 

La figure 10 presente une comparaison entre des prof tls de vitesse mesures et les profils predits par ie 
modele decrit ci-dessus. Les profils representes correspondent aux conditions suivantes : debit total de 4000 
40 barils/jour (26,5 m 3 /h), et inclinaison 9 de respective ment 5°, 25° et 45°. On constate une concordance conve- 
nable entre les valeurs experimentales et les resultats obtenus a partir du modele. 

Outre le modele de vitesse pour les valeurs de P importantes (expression (11)), un modele de vitesse a 
ete developpe pour P petit (P<1). Dans ce cas, les gradients d'ecoulement sont petits et uniformes. Le modele 
est def ini comme suit: 

*5 v(z 2 ) = V m + 6V* 2 (Z2 - 1/12) + 2V* 2 .Z - 2V* 2 .e2P<z- (15) 

ou V* 2 est une constante egale a 0,3 m/s environ. 

A partir des modeles specif iques des expressions (11) et (15), le modele general suivant a ete etabli: 

v(z 3 ) = e- ww.v(z t ) + (1 - e - p/10 ).v(z 2 ) (16) 
Pour P grand, I'expression (16) se ramene a I'expression (15). Pour P petit, I'expression (16) se ramene 
50 a I'expression (11). 



Revendications 

1. Procede d'analyse de I'ecoulement dans un puits incline, recoupment comprenant du petrole et de i'eau, 
caracterise par le fait que, pour une profondeur ou Ton dispose d'au moins une mesure d'un parametre 
d'ecoulement, on etablit a partir de cette mesure, connaissant i'inclinaison 9 et le diametre D du puits a 



7 



# 



EP 0 519 809 A1 



la dite profondeur, un prof il de variation de ce parametre le long d'un axe z'z def in i par le diametre de la 
section transversale du puits situe dans un plan vertical, en utilisant un mod el e prevoyant une couche de 
melange eau-petrole turbulente. 

Procede selon la revendication 1, dans lequel le dit parametre est la proportion statique d'eau a et on 
etablit le prof il de variation ct(z) en ajustant le modele a partir d'au moins une mesure locale du dit para- 
metre. 

Procede selon la revendication 2, dans lequel on determine a partir du prof il ct(z) la valeur moyenne H w 
de la proportion statique d'eau sur la section transversale consideree. 

Procede selon la revendication 1 , dans lequel on etablit le prof il de variation en ajustant le modele a partir 
d'une mesure de la valeur moyenne du dit parametre sur la section consideree. 

Procede selon Tune des revendicattons 1 a 4, dans lequel la couche de melange turbulente a une dimen- 
sion X selon I'axe z'z donnee par la relation 

KsinO 

dans laquelle K est une constants pour un diametre D donne. 

Procede selon la revendication 5, dans lequel le modele prevoit un gradient constant pour la variation de 
la proportion statique d'eau a dans la couche de melange, selon la relation 

ou Zc est la position sur I'axe z'z du centre de la couche de melange. 

Procede selon Tune des revendications 1 a 4, dans lequel le modele prevoit que la proportion statique 
d'eau a varie suivant I'axe z'z conformement a la relation : 

a(Z) = 1 - ®2 pz Pour z < minfe, D/2} 

zsnp 

oCi 

-Z = z/D 

- X est une constante pour un diametre D donne 

■ p = o^ XS ' n9 v ou a m est la valeur moyenne de a suivant le prof il, 
2(1 - 

- Ze est def inie par <x(Ze) = 0 

- min{z e , D/2} vaut D/2, ou Zq si z« est plus petit que D/2, 

et dans le cas ou z e est plus petit que D/2, le prof i I comprend en outre une portion def inie par : 

a = 0 pour ^ < z < D/2. 

Procede selon la revendication 1, dans lequel le dit parametre est la vitesse, et on etablit un prof il de va- 
riation de la vitesse v(z) en ajustant un modele a partir d'une mesure de la valeur moyenne de la vitesse 
sur la section consideree. 



8 





9 



* 



EP 0 519 809 A1 



FIG. 3 



-as 



B4 




B2 



S3 



FIG. 4 



N 




EP 0 519 809 A1 



FIG. 5 



EMT#E£ D£S 
VALEt/KS ZX£ 
I), &, K 



54 



£MT&££ D£S 
YALE///? J?£ M£SU#£ 
2> t , <*l , St 



52 



2?£S Ct?OXI>OW££S 
2?£S SO/VD£S 



/H/T//JL/ SAT/OS/ J>£ 

I0y£/W£ 23£S o(£ 



M/A//M/SA T/OW J>£ 



//VT£G*/lT/OA/ Z>(/ 



53 



56 



11 



* 



EP 0 519 809 A1 



i 



FIG. 6 



£AfT/?££ Z>£S VAL£Uf?S 



/N/TMLIS4T/0A/ J?£ of^ 



sor la sscr/o/y 



ou/ 



S2 



€3 




65 

/ 



MOZ)/F/£& 

Z3 



Af0A/ 



12 



# 



EP 0 519 809 A1 




FIG. 7 
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FIG. 9 
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